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Abstract 

We investigate universal features of the jamming transition in granular materials, colloids 
and glasses. We show that the jamming transition in these systems has common features: 
slowing of response to external perturbation, and the onset of structural heterogeneities. 

1 Introduction 

Many studies in physics concern the onset of instabilities into turbulence or some form of 
chaotic break up, the simplest being laminar flow into turbulent flow, but the collapse of 
civil engineering structures is much studied. There are not many studies of the reverse 
phenomenon, where some chaotic system becomes more regular, or even stops altogether. 
The jamming phenomenon exists in quantum as well as classical systems, but in this article 
we will discuss only the latter, and indeed systems which are readily accessible to intuitive 
understanding. 

The approach of this article will be to list the most common jamming systems, and 
see if there exist simple physical or mathematical systems which lead to interpretable 
experiments and soluble theories. In this paper, we consider mostly our own work, or 
studies that we know of in detail. The word "jamming" derives from systems coming to 
rest, and an archetypical problem is in the flow of granular or colloidal systems. We will 
try to answer the following questions: 

• What is the jammed state? 

• How does the jamming transition occur? 

• What are the common features of the jamming transition in granular materials, 
colloidal suspensions and glasses? 

In this paper, we confine ourselves to the discussion of the jamming transition in 
granular materials, colloids and glasses. We will not comment on other systems that 
exhibit the jamming transition, which are many: foams, vortices in superconductors, field 
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lines in turbulent plasma etc. All these systems throw light on one another and exhibit 
universal behaviour. We will demonstrate that the jamming transition in these systems 
can be characterized by the slowing of response to external perturbations, and the onset 
of structural heterogeneities. 



2 Granular Media 

The jamming transition in dry granular media is a very common phenomenon which 
can be observed in everyday life. Stirring jammed sugar with a spoon can be hard but 
manageable. A jam within a silo can cause its failure. Nonetheless, such an ubiquity does 
not make the problem less difficult. In this section we will not even attempt to analyze the 
dynamical problem of a granular flow coming to rest, and the onset of jamming. Instead, 
we will try to characterize the static jammed state of a granular material. The difficulty of 
a theoretical analysis lies in the fact that dense granular media can expand upon shearing, 



a well known Reynolds dilatancy, that can be viewed as a counterpart of jamming [15, 25]. 
Consider a cylindrical vertical pipe in the following cases (see Figure [j]) : 

• Rough walls, filled with pieces of twisted wire. This wire will entangle and not flow 
at all, indeed one does not need the pipe. 

• Rough walls, with rough, approximately spherical particles. These can flow under 
certain conditions, but can also jam; it is a well studied problem in chemical engi- 
neering. 

• Smooth walls, smooth spheres. This is difficult, perhaps impossible, to jam. 




Figure 1: Types of granular materials 



If one asks for the simplest granular material in the above list, it must be the third one 
if we model the granular material as an assembly of discrete rigid particles whose interac- 
tions with their neighbours are localized at point contacts. Therefore, the description of 
the network of intergranular contacts is essential for the understanding of the stress trans- 
mission and the onset of the jammed state in granular assemblies. The random geometry 
of a static granular packing can be visualized as a network of intergranular contacts. For 
any aggregate of rigid particles the transmission of stress from one point to another can 
occur only via the intergranular contacts. Therefore, the contact network determines the 
network of intergranular forces. In general, the contact network can have a coordination 
number varying within the system and different for every particular packing. It follows 



then that the network of intergranular forces is indeterminate i.e. the number of unknown 
forces is larger than the number of Newton's equations of mechanical equilibrium. Thus, 
in order for the network of intergranular forces to be perfectly defined, the number of equa- 
tions must equal the number of unknowns. This can be achieved by choosing the contact 
network with a certain fixed coordination number. In this case the system of equations 
for intergranular forces has a unique solution and the complete system of equations for 
the stress tensor can be derived. This is the simplest model of a granular material. If 
this cannot be solved, one will be left with empiricism. This is from the point of view 
of physics; it is not a good philosophy for engineers. The geometric specification of our 
system is as follows: we will need z contact point vectors TZ a ^, centroid of a grain a, 
R a , f a/3 the vector from R a to TZ a ^, and the distance between grains a and /?, R a @ (see 
Figure ^). 




Figure 2: Detail of two grain contact 



Grain a exerts a force on grain j3 at a point TZ a ^ = R a + r a/3 . The contact is a point 
in a plane whose normal is n al3 . The vector R a is defined by: 

R a = ^—, (1) 

so that R a is the centroid of contacts, and hence the relation: Y^a"^ = 0- We note 
that z is the number of contacts per grain, and ^2 a means summation over the nearest 
neighbours. We define the distance between grains a and (3 

R^ = r Pa - r afi , (2) 

Hence R a , r a/3 and n a ^ are geometrical properties of the aggregate under consideration 
and the other shape specifications do not enter the equations. In a static array, Newton's 
equations of intergranular force and torque balance are satisfied. Balance of force around 
the grain requires 

£/^ = <?f, (3) 

P 

fi* + ff a = 0, (4) 



where g a is the external force acting on grain a. 
The equation of torque balance is 

Y,^iffr^ = Ct. (5) 

P 

Friction is assumed to be infinite [17]. It can be verified that, for the intergranular forces 



in the static array to be determined by these equations, the coordination number z = 3 
in 2-D and z = 4 in 3-D is required. The microscopic version of stress analysis is to 
determine all of the intergranular forces, given the applied force and torque loadings on 
each grain, and geometric specification of a granular array. The number of unknowns 
per grain is zd/2. The required force and torque equations give d + constraints. 
The system of equations for the intergranular forces is complete when the coordination 
number is z m = d + 1. In addition, the configuration of contacts and normals is only 
acceptable if all the forces are compressive. If tensile forces are allowed, then we would be 
studying a sandstone (an interesting problem, but not a subject for a study of jamming). 
Many investigators do not believe that z will be this small, and would invoke the kind 
of arguments used in bridgework, where the same problem arises when extra spares lead 
to too few equations for solution. However the simplest assumption is to assume that if 
the geometry gives more than contacts, some will contain no force. At all events we 
can restrict ourselves to systems where there really are ^ contacts per grain The main 
question is: can one observe the jammed state in this simple model where the static stress 
state can be determined? We define the tensorial force moment: 

s£ = £ifV> ( 6 ) 

p 

which is the microscopic analogue of the stress tensor. With Cf = 0, will be symmetric. 
To obtain the macroscopic stress tensor from the tensorial force moment, we coarse-grain, 
i.e. average it over an ensemble of configurations: 

N 

°*<f) = <£ S°5(r-R a )). (7) 

a=l 

The number of equations required equals the number of independent components of a 
symmetric stress tensor aij = aji, and is rf ^ 2 + ^ . At the same time, the number of 
equations available is d. These are vector equations of the stress equilibrium = 

which have their origin in Newton's second law. Therefore we have to find di ~ d ~ 1 ^ equations, 
which possess the information from Newton's third law, to complete and solve the system 
of equations which governs the transmission of stress in a granular array. 

Given the set of equations (|| — ||) we can write the probability functional for the 
integranular force as 



p{f i ap } = mj2fi aP -9?) 

P 

x SU? + ft) , (8) 



where the normalization, J\f, which is a function of a configuration, is 



A/ 1 



(9) 



The probability of finding the tensorial force moment on grain a is 



a,/3 (3 



a 3 



(10) 



where J Vf a @ implies integration over all functions f a @, since all the constraints on f a @ 
have been experienced. We assume that the z = d + 1 condition means that the integral 
exists. 



It has been shown |18|, [I9j, |20[ that the fundamental equations of stress equilibrium 
take the form 



VjCTjj + V jVfcV m KijklOlm + ... = g 



(11) 



PijkPjk + V 'jTijkl<Jkl + VjVlUijkl<Tkm + ... = . (12) 

However, there are difficulties with the averaging procedure [18] which is highly non-trivial. 
The simplest mean-field approximation gives the equation an = o<i2 for the case of an 
isotropic and homogeneous disordered array [O] • Though this equation gives the diagonal 



stress tensor aij = pdij, it is not rotationally invariant. The only linear, algebraic and 
rotationally invariant equation is Trcr^ = 0. However, this equation cannot be accepted, 
for a stable granular aggregate does not support tensile stresses. We believe that, at least 
in the simplest case, the fundamental equation for the microscopic stress tensor should 
be linear and algebraic (because of the linearity of Newtons' second law for intergranular 
forces). In this paper we offer an alternative way which is considered below. 

The leading terms of the system of equations (111, 12) arise from the system of discrete 



linear equations for S?f |18|, [19] 



SijMjlR k 



a/3 



(13) 
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2S& tan 9 a 



(14) 



and if tan a has an average value tan < 



en 



0"22 



2o"i2 tan ( 



(15) 



which is known as the Fixed Principal Axes equation [21, 22, 23 



], and has been used 

with notable effect to solve the problem of the stress distribution in sandpiles. For a 
homogeneous and isotropic system, the averaging process gives the stress tensor = p 5ij 
which is simply hydrostatic pressure, as is to be expected. Rotating Sfj by some arbitrary 
angle 9 one can easily that (01) constraints the off-diagonal components of Sfj to be zero. 
The system of equations (|l3jl4]) is solved by Fourier transformation and the macroscopic 
stress tensor is obtained by averaging over the angle 9 



ia n (k) = (Sn(k)) e 



+ 92 (k% 



(16) 



• /A /Q ,A\ g 2 (fc| + 3fcffc 2 ) + gl(M| ~ fcj) / 17 x 
W22(«) = {S22{k))e = (17) 

| A; | 

■ /A /A \ (gife -g 2 ki)(k1 - k\) 

ia 12 (k) = {b 12 {k)}e = (18) 

\k\ 

where \k\ 2 = k\ + h\ and <Jij(f) = / <Jij(k) e lkr d?k. By doing the inverse Fourier 
transformation one can see that the macroscopic stress tensor is diagonal. There must 
also be constraints on the permitted configurations (due to the absence of tensile forces) 
which are not so easily expressed, for they affect each grain in the form 

SgM^Rfuf > (19) 

which has not yet been put into continuum equations other than Det a > and Tr a > 0. 
However, this condition must be crucial, for without it jamming becomes a very ubiquitous 
phenomenon provided the density of grains is such that they are all in contact. This 
example shows the utility of simple models. Apply this type of grain in a vertical pipe 
(the grain is a thick tile, see Figure |||). 




Figure 3: Jammed pipe 



This system is clearly jammed, but the problem is how it can get into this configuration. 
Another limiting case is to think of a sphere with many spines pointing in the radial 
direction. These spines mean that spheres can only approach or retreat along the direction 
joining their centres, and the only other motion permitted to a group of two or more is 
that of rigid rotation of the group. As soon as a line (in 2-D) or shell (3-D) of these 
objects occurs, they jam; so they always jam. Again, although this is a possible material, 
most materials fall into the classes above. So there are trivial jamming problems, but it 
is proving difficult to produce an effective analytic theory for intermediate materials. It 
is natural then to use computer simulations, and there is a notable literature in existence 
|25|| . However, rather than comment on this literature for granular materials we move to 
the related problem of colloidal flow which offers a natural route from grains to glasses. 



3 Colloidal Suspensions 



The simplest model of a colloid which exhibits the jamming transition is the sheared sus- 
pension of hard monodisperse spheres, interacting hydrodynamically through a Newtonian 



solvent of viscosity ijq. Such a system at equilibrium is characterized by its volume fraction 



<f> v . The behaviour of this system at equilibrium is well known 1 33 1 . With increasing (f> v , the 
system crystallizes, with phase coexistence occurring between <p v ~ 0.50 and <p v ~ 0.55. 
If simple shear is applied, we need one other parameter, which is the Peclet number 

Pe = (20) 
k B T 

where d is the particle diameter, 7 the shear rate and T the temperature. The Peclet 
number gives the ratio of shear forces to Brownian forces. Pairwise hydrodynamic inter- 
actions between the neighbouring particles can be divided into squeeze terms along the 
line of centres and terms arising from relative shear motion. To leading order the squeeze 
hydrodynamic force on a particle is given by the well-known Reynolds lubrication formula 

f* = - E^{("i - + °( ln ^) (21) 

where the sum is over nearest-neighbour particles j, hij is the gap between the neoghbour 
surfaces with the unit of distance the particle diameter, riij is the unit vector along the 
line of centres i to j and Vi, Vj are the particle centre velocities. 

In the absence of other interactions, Brownian forces are left to control the approach 
of particles. If conservative (steric or charge) forces are present, these control the gaps 
between particles, and dominate over the Brownian forces at high shear rates. When one 
studies the shear flow of such a system, it exhibits thickening effects: a rise of viscosity with 
increasing shear rate (for a review, see |34[]). At volume fractions approaching random close 
packing, (j)RCP ~ 0.64, discontinuous thickening with a large jump in viscosity occurs at 
a critical Pe. However, at a lower cp v and lower Pe a more continuous rise can be observed 
(see the Figure |] below). Hence, the jamming transition in this simple system can be 
either continuous or discontinuous, which depends sensitively on the volume fraction. 

Jamming of hard sphere colloids 

Simulations ol 700 particles at different volume fractions, ®, 



Figure 4: Relative viscosity of hard sphere colloids versus time, for various volume fractions 



This happens in various experimental systems but in particular, in those whose parti- 
cles, by polymer coating or surface charges, do not flocculate (if they want to stick together 
then jamming is not an obvious phenomenon). Colloids with repulsive interactions exhibit 
shear thinning i.e. the viscosity decreases as the shear rate increases. The presence of ag- 
gregating forces can greatly alter the shear thinning. We will discuss the regime of shear 
thickening (which we call later the jamming transition) in the simplest model, gradually 
incorporating various interactions. The physical picture can be obtained by combining 

]. The simulation technique for 



theory and computer simulations [35, ^g, 37, £38], |3^, 40 



particles under quasi-static motion determined by a balance of conservative and dissipative 
forces has been proposed in [35|. The motion of N colloidal particles, immersed within a 
hydrodynamic medium, is governed by an equation of quasi-static force balance 



- 1Z(X) ■ V + F C (X) + F B (t) = (22) 

where X represents the 67V particle position coordinates and orientations, 1Z is a 6N x 6N 
resistance matrix and V = dX / dt is the particle velocity. The terms Fc and Fb represent 
conservative and Brownian forces. The effects of inertia are ignored i.e. the Reynolds 
number is small 

Re = »1* « 1 (23) 

Bearing in mind that for a typical colloid ^ ~ 10 7 , it follows that it is possible to achieve 
Re < 1 and Pe > 1 simultaneously. Equation (f22|) can be solved on a computer p|, 13(1. 



Although simulation shows exactly what is happening in the sense one can see every sphere 
flowing its path (see Figures |],||), it is still not agreed as to whether the phenomenon of 
thickening is an order-disorder transition p2|, |43| or whether it is due to the development 



of clusters of particles along the compression axis [44 



Figure 5: Snapshot from a 3-D simulation of 2000 hard spheres at 4> v = 0.56. All particles 
are shown. 



A microscopic kinetic theory for the origin of the increased bulk viscosity at high (f> v and 



Pe = oo has been recently proposed [38, |40[| . It attributes the viscous enhancement to the 
presence of hydrodynamic clustering, i.e. incompressible groups of particles which lie near 
the compression axis. The rigidity of clusters is provided by divergent lubrication drag 
coefficients. The theory predicts a critical cp v above which jamming occurs. This model 
gives a flow-jam phase transition at any strain and may be of more general applicability 
(e.g. colloids with conservative and Brownian forces). Assuming the cluster length to be 
additive on collision, one can obtain the standard Smoluchowski aggregation equation 



dXk 
dt 2 



^ oo 

2 X! Xi Xj — X i+ j \Si + j t k — 5i,k — Sj t kj (24) 



This equation relies upon a mean field approximation, i.e. it is assumed that each cluster 
is embedded in an average flow, composed of the other clusters, "spectator" particles and 
solvent. Eqn (p3) governs the evolution of the concentrations X& of clusters of size k 




Figure 6: Snapshot from a 3-D simulation of 2000 hard spheres at <j) v = 0.56. Only 
particles who have separations h < 10 d are shown. The flow direction is left to right. 



monomers as a function of time, given that clusters can collide, aggregate and break up 
at rates dependent on their sizes and specified by the aggregation kernel K^j. The math- 



ematical procedure of solving Equation (24) has been reported in 4C]. We discuss 



the generic structure of the theory. The rheology of this system results from a compe- 
tition between a binary aggregation process and single cluster breakup. The paradigm 
of hydrodynamic clustering modelled in terms of aggregation-breakup laws provides the 
qualitative physical picture of the jamming transition. At low <f> v breakup dominates. The 
system achieves a steady state with a population of clusters whose average size increases 
and a viscosity which increases as the volume fraction rises. In the vicinity of some volume 
fraction cp c , average cluster size diverges, and so does the viscosity 



7/ oc % 



C 1 - x)~ 2 < 25 > 



Above (f> c the flow is transient and dominated by aggregation. The average cluster size 
diverges after a certain strain, which falls as <p v is further increased. The system undergoes 
the jamming transition and can never reach a steady state. This model predicts C ~ 0.49, 
although experiments suggest a divergence of bulk viscosity at a higher volume fraction 
(which is still less than 4>rcp)- This difference is due to the presence of conservative 
interactions in the real hard core colloid. In conclusion, we briefly discuss the jamming 
behaviour of long flexible chain polymers. A good way to understand this is to plot viscosity 
against the molecular weight of the polymer in the molten state. The molecular weight is 
equivalent to the length of the chain. The short chains slide past each other in many ways, 
giving a linear relation between length and shear viscosity, but an entanglement crisis 
occurs at a critical molecular weight M c , when the polymer can only wriggle in Brownian 
motion up and down and "out" of a tube formed by its neighbours in reptation. Whether 
there is sharp transition is not established in a difficult slow experiment, but all agree 
that the dependence of the viscosity jumps to M 3A , a colossal change which at first sight 
is just like a jamming, but very slowly the melt flows flows due to reptation. Polymers 
very easily form glasses and it only takes a fall in temperature to make the melt solidify, 
although the specific heat shows that although reptation (a slow translatory motion) has 
ceased, there is plenty of other movement taking place, but movement which averages over 
time to zero. A further fall in temperature destroys this motion, and one enters a state 
like a conventional glass. 



4 Glasses 



Studies of glasses seem to fall into two camps. Those taking continuum field representation 
of the material and solving under statistical thermodynamics conditions, principally mode- 
coupling methods, and alternatively the translation, into comparatively simple equations 
and simple physical models. The former has given intuitive thinkers a hard time because 
it is not obvious in those cases, such as the polymeric glasses discussed below, where one 
should have a picture of what goes, what indeed is going on. The modelling approach is 
weak; for example, in an assembly of packed spheres where an enormous number of motions 
is possible, but strong in polymer glasses where the motion is obvious and the jamming 



of motion is the cessation of the centre of mass diffusion. The mode-coupling theory[46] 
of the glass transition for simple liquids has an ergodicity to non-ergodicity transition. 
It starts with an equation of motion for the density correlation function, which is an 
integro-differential Langevin equation, with a non-linear memory kernel. This non-linear 
memory term governs the transition to the non-ergodic state, and can be characterized by 
measurements of the dynamic structure factor. We refer the reader to the literature f46[ l for 
details of this approach. The glass transition is not limited to special types of materials. 
Every class of material can be transformed in an amorphous solid if the experimental 
parameters are adjusted to the dynamics of the system. Consider therefore, two extreme 
examples. Consider the system consisting of spherical molecules such as rare gases. The 
hopping time of the spheres is very short, and the dynamics extremely fast. Nevertheless, 
it has been shown that such fluids undergo a glass transition from the super cooled melt, 
if one cools the system with a quenching rate of q oc 10 12 Ks _1 , when the molecules jam 
effectively to rest. On the other hand, consider polystyrene which consists of very large 
molecules. The dynamics of such a system is much more complicated in comparison to 
spheres because there are many degrees of freedom. It is most significant that the centre 
of mass diffusion of a single molecule is small. We want to stress characteristics of the 
glass transition in general. The most significant points we want to discuss are: 

• The divergence of the transport or inverse transport properties, such as viscosity, 
inverse diffusion coefficient, and relaxation times. 

• The extreme broad relaxation phenomena of the stress, modulus etc. 

• The quantitative definition of the term cooperativity. 

• Influence of external parameters on T g . 

It has been recognised that the relaxation time follows the Vogel - Fulcher Law (VF): 

r oc (26) 

This law has many names. In polymer physics it is often called the Williams - Landel 
- Ferry ( WLF) law jl?]] . The law is not valid over the whole temperature range (see Figure 
0). The divergence given by ( p6[ ) comes at To, where To is a temperature below the freezing 
temperature T g , and the empirical rule is given by: 

c 

T ~ T g - (20 ~ 30)°. (27) 

The physical meaning of T g is still unclear, and we want to try to clarify this point. The 
VF law is much stronger than the critical slowing down in phase transition phenomena 
where, by scaling arguments, the relaxation time is given by 




Figure 7: Vogel-Fulcher Law 



r ~| T - T \~ vz (28) 

where v is the correlation length exponent i.e. £ ~ (T — Tq)~ u and z the dynamical 
exponent 

r ~ i z . (29) 

There have been attempts to fit data for freezing transitions (glass transition, spin glass 
transition) by Eqn.(^), and it turned out that vz is very high vz ~ 10, ...,20 which 
seems very unnatural. This again indicates a physical significance to the VF law. Another 
peculiar point lies in the relaxation properties. An empirical law was found long ago by 
Kohlrausch and in the early seventies recovered by Williams and Watts in their studies of 
the broadening of relaxation processes [f5ij |. For example, measurement of the dielectric 
constant shows a much larger half width in the imaginary part compared to the Debye 
process. Empirically this is described by the Kohlrausch - Williams - Watts (KWW) law 
(see Fig D 

e" (co) 




Figure 8: Kohlrausch- Williams- Watts Law 



where 4>{t) can be any quantity which relaxes, i.e. 



(30) 



_ e(t) - 6(00) 
W " e(Q) _ e(oo) W 

However, there has been no convincing explanation for a unique value for (3. We doubt 
that there is more to the physics of (^) than a wide relaxation spectrum due to different 
physical processes, and a more relevant question is how both the KWW and VF laws link 
together in general, and what is the relationship to cooperativity. It is believed that the 
glass transition exhibits a large amount of cooperative motion as the system is close to T g . 
This might be indicated as well by the VF law, which is the crossover to the Arrhenius 
behaviour right at T g but with an extremely high activation energy (see Figure This 
activation energy is so high that it could hardly be attributed to only one molecule, and 
the phenomenological interpretation is that there are cooperative regions of some linear 
size diverging at some temperature 

£ ~ (T - T g ) power . (32) 

Another quite general question is the state of time-temperature superposition principle. 
This says that if one measures a physical quantity, D (t), at some temperature T, and if the 
measurement is repeated at some temperature Ti, the quantity D(t) can be resolved by a 
"shift factor" ax, i.e. D(T,t) is not a function of two variables but only of a combination 
of both: 



D(T,t) = D(— ) (33) 

where ax is often given by 



near T g , and 



aT = T-T g + C 2 (34) 



SE 5E 

a T = ~Y - ^7 ( 35 ) 

far from T g . Hence ( |34| ) is of the Volger-Fulcher form. 
Polymeric glasses offer two challenges: 

• There is a clear intuitive picture of what is happening. The "tube" closes at its ends, 
or contracts at "entanglement" points. 

• The molecular weight offers, as with viscosity, a new degree of freedom, and hence 
new laws emerge. Any theory of glass must encompass VF, KWW, and the experi- 
ments we now describe. 

For polymers, the glass transition temperature depends on the molecular weight, and an 
empirical rule is given by Flory and Fox: 



T g (L) = T g (oo)-^^, (36) 



constant 
T 

where L is the length of the molecules i.e. the molecular weight. The agreement of ( |36| ) 
with experiments is not particularly good, but it gives an estimate of T g (L). Hence ( |36| ) 
tells us that there is not a very significant dependence on the molecular weight, unless L 
is small. This point should be investigated using the knowledge of polymer dynamics in 



melts which has recently emerged. Another empirical law we want to discuss is the mixing 
rule in plastification. Mixing two glass forming polymers together, the new T g is often 
given by 



' = ZL + ZL (37) 
T 9 (mix) T gi T g2 

to zeroth order. The (pi is the volume fraction of the i'th polymer. We have discussed so 
far the most important experimental results. Clearly there is a need for a solvable model 
in the framework of which the VF and KWW laws can be derived. Any fundamental 
theory of glass transition should relate the mobility of a molecule within the cage to 
mobilities of molecules forming the cage. These mobilities of surrounding molecules are 
coupled with those of their neighbours etc., and therefore in general there is no small 
parameter that can justify the decoupling approximation. The small molecule model of a 
glass transition involves cages, but polymers are simpler because they have a tube, and 
a tube of such a model can be modelled as a straight tube, indeed restrictions along the 
tube can be included, but are not here for simplicity. The reason why we take rods was 
the advantage of simple geometry, no internal degrees of freedom, and very slow dynamics, 
so that we have no problems with high quenching rates, so we do not have to worry about 
thermodynamics. It is obvious that the small width to the length ratio and the sufficiently 
large number of neighbouring rods justify the decoupling procedure. If the solution of the 
rods is dense, the concentration c < d 2 L {d is the diameter of the rod, L its length and c 
the concentration) and severe constraints are acting on the rods. For example, they cannot 
move rotationally, and they can only make progress along their length. Such a solution 
has been called entangled. Suppose now we have such a solution of highly entangled rods. 
A rod can slide between the entangling rods until it meets rods which block it (see Figure 
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Figure 9: "Jammed" rod 



The motion of a rod will then be like a particle diffusing along a line but meeting gates 
which open and close randomly through thermal fluctuations. If no barriers we present 
the probability P(x, t) of finding the test rod at x (which is the coordinate of the rod down 
the tube) and time t satisfies the simple diffusion equation 




This has the solution 



/oo 
G Q {x,x';t,t')P {x',t')dx' (39) 
-oo 

where Pq{x' ,t') is the initial probability function and Go(x,x';t,t') is the standard Green 
function of the diffusion equation. Suppose now , a reflecting barrier is placed at position 
R at time tR and removed at some time tQ. Then using the method of images to find the 
Green function it is straightforward to calculate P(x,t) [|49| . The same method can be 
applied when there are many barriers appearing and disappearing along the path of the 
rod. After neglecting correlations between barriers one can see that the rod (1) can only 
diffuse if 



D = D (l - a) (40) 



where a ~ e(cDL 2 )2, cDL 2 being the Onsager number. Eqn.(|40|) gives D = if 
a = 1 i.e. ecDL 2 = 1. This is the jammed state. The result can be modified to 



include cooperativity. The complete solution has been obtained in [49, |5(| by mapping the 
problem onto the self-avoiding walk problem. The VF law appears when summed over n 
(where n is the number of rods that loops are made out of) 



D ~ D exp ( - a% ) (41) 

\ 1 — Oil' 



where the parameter a\ contains generalised constants and a minimum loop size, and 
a 2 = ol — Oil. It can be shown that the number of rods moving cooperatively in the loops 
are given by 

n ~ (1 - a 2 y 2 (42) 

and the size of the loop is therefore 

e - nL 2 ~ ^ (43) 

(i - u 2 y 

Thus £ diverges if a 2 — > 1, as the phenomenological interpretation requires. So far we 
have used a to denote the expression cdL 2 , a combination of constants well known in 
liquid crystal theory since the work of Onsager. However, its role in our theory is much 
more general; for example, d can be temperature dependent through the fact that Van der 
Waals forces appear in the form e~^r, whereas hardcore forces do not contain T. Thus, 
at the level of a model one can regard a = 1 at T = T g . Until one studies a detailed 
physical case, there is no purpose in false verisimilitude. Thus the above model now says 
that the VF law is a direct consequence of cooperativity. Turning now to the relaxation 
behaviour, it can be shown that for the rod model, the stress relaxation follows, to first 
order, the law 

a(t) ~ e~ ( ^ (44) 

Cooperativity as indicated by (^l]) does not change ( f44| ) drastically, giving rise to loga- 
rithms 

a{t) ~ e~ ( " )!? [I + const, log -J (45) 

This leads us to the conclusion that cooperativity is responsible for the VF behaviour, 
but not for the relaxation phenomena. Further study of this model allows a derivation 



of (36,37). A final interesting feature of the tube model is that the tube itself must be a 
random walk of step length a (in the rheology literature this is called the primitive path) 
and the "rod" of the preceding discussion is of length a for a very long polymer, and L for 
a long, but not very long molecule. The freezing of large scale motion is the freezing of 
motion on the scale of the primitive path step length, and the long term reptative motion 
of the whole chain is not a vital constitutive of the glass temperature. Hence T g is a 
function of a and the density of the material, and also additional parameters (e.g. chain 



stiffness), i.e. T g 



T g (a{c, L), 



But a is a function of the density itself, so that the 



effect of a diluent acts on both a and c, and we expect a different concentration dependence 
of T g above M c . T g dependence on L is roughly given by Figure 1C. 



T (a, c) 




* crossover 



Figure 10: T„ as a function of L 



This simple model has been modified for arrays of rigid rods with fixed centres of 
rotation [51]. The rods attached to the sites of the cubic lattice can rotate freely but 



cannot cross each other. It is important to note that the glass transition in this system 
is decoupled from the structural transition (nematic ordering) and the only parameter is 
the ratio of the length of the rods to the distance between the centers of rotation. The 
Monte-Carlo study in 2-D shows that with increasing the parameter of the model a sharp 
crossover to infinite relaxation times can be observed. In 3-D the simulation gives a real 
transition to a completely frozen state at some critical length L c . 

4.1 An Important Analogy 

Recent crucial experiments on granular materials 52, [53| show that external vibrations 



lead to a slow approach of the packing density to a final steady-state value. Depending on 
the initial conditions and the magnitude of the vibration acceleration, the system can either 
reversibly move between steady-state densities, or can become irreversibly trapped into 
metastable states; that is, the rate of compaction and the final density depend sensitively 



on the history of vibration intensities that the system experiences (see Figure 11). 

The function which has been found to fit the ensemble averaged density p{t) better 
than other functional forms, is J53|] : 

Pit) =Pf ^7 FT ( 4 6) 

HKJ rt 1 + B log(l + 1) V ; 

where the parameters pf, A pj and r depend only on T. This is an analogue of a glass 
forming material where the lower curve corresponds to the situation where the quenching 
speed inflicted on the glass is faster than the speed at which the glass relaxes back to 



P A 
P(X=0) - 



P(t=0) 




r 

Figure 11: Dependence of the steady-state packing density on the tapping history (Nowak 
et al.). Experimental values of density packing fraction are in the following correspondence 
with model parameters: p(X = 0) = — 0.64, p(t = 0) = po = — « 0.58 and p(X = 
00 ) = [vq+vi) ~ 0.62. The vibration intensity is parametrized by V = | 



equilibrium. The special feature of the granular material is that it has no power of its 
own to proceed to the equilibrium state, so that every aspect can be studied without 
having to worry about the fact that the true glass is always seeking equilibrium. It is 
worth considering a simple model for the density response to external vibrations |{54| . If 
we assume that all configurations of a given volume are equally probable, it is possible 



to develop the formalism |5q, Kq, 57, 58, 5£, pOl analogous to conventional statistical 



mechanics. We introduce the volume function W (the analogue of a Hamiltonian) which 
depends on the coordinates of the grains, and their orientations. Averaging over all possible 
configurations of the grains in real space gives us a configurational statistical ensemble, 
which describes the random packing of grains. An analog of the microcanonical probability 
distribution is: 



s 



p = e -*5(V - W) . (47) 
We can define the analogue of temperature as: 

dV 

x = w (48) 

This fundamental parameter is called compactivity |^]. It characterizes the packing of a 
granular material, and may be interpreted as being characteristic of the number of ways it 
is possible to arrange the grains within the system into a volume A V, such that the disorder 
is AS. Consequently, the two limits of X are and oo, corresponding to the most and 
least compact stable arrangements. This is clearly a valid parameter for sufficiently dense 
powders, since one can in principle calculate the configurational entropy of an arrangement 
of grains, and therefore derive the compactivity from the basic definition. We will use the 
canonical probability distribution 

Y-W 

p = e^r, (49) 

where A is a constant which gives the entropy the dimension of volume. We call Y the 
effective volume; it is the analogue of the free energy: 



Y 
' XX 



/W(a) 
e~— d(all), (50) 

dY 

v = y-x-. (5i) 



Examples of volume functions for particular systems can be found elsewhere [57, 58, 
We consider the rigid grains powder dominated by friction deposited in a container which 
will be shaken or tapped (in order to consider the simplest case, we ignore other possible 
interactions, e.g. cohesion, and do not distinguish between the grain-grain interactions in 
the bulk and those on the boundaries). We assume that most of the particles in the bulk 
do not acquire any non ephemeral kinetic energy, i.e. the change of a certain configuration 
occurs due to continuous and cooperative rearrangement of a free volume between the 
neighbouring grains. The simplest volume function is 

W = v + (ui - vq){hI + (52) 

where two degrees of freedom [i\ and fi2 define the "interaction" of a grain with its nearest 
neighbours. If we assume that all grains in the bulk experience the external vibration as 
a random force, with zero correlation time, then the process of compaction can be seen as 
the Ornstein-Uhlenbeck process for the degrees of freedom fa, i = 1,2 |pl| . Therefore we 
write the Langevin equation: 

dm IdW r^ fM 

-JT + = vDfi(t) (53) 

dt v djii 

where {fi(t)fj(t')) = 28ij5(t — t'), and v characterizes the frictional resistance imposed on 
the grain by its nearest neighbours. The term fi(t) on the RHS of ( |53"1) represents the 
random force generated by a tap. The derivation of this gives the analogue of the Einstein 
relation that v = (XX)/ D. If we identify / with the amplitude of the force a used in the 
tapping, the natural way to make this dimensionless is to write the "diffusion" coefficient 
as : 

D=(-) , 54 

that is we have a simplest guess for a fluctuation-dissipation relation: 

'a\2 i> 2 uj 2 



«=rr^ (55) 



-9' v 

where v is the volume of a grain, u the frequency of tapping, and g the gravitational 



acceleration. The standard treatment of the Langevin equation (|53|) is to use it to derive 
the Fokker-Planck equation: 

dP / d 2 d \ 

where Dij = DSij and jij = j5{j. This equation can be solved explicitly |34|, [6l[| . We 
can calculate the volume of the system as a function of time and compactivity, V(X,t). 
Though this is a simple model, it is too crude to give a quantitative agreement with 
experimental data; however, it gives a clear physical picture of what is happening. It 
is possible to imagine an initial state where all the grains are improbably placed, i.e. 
where each grain has its maximum volume v\. So if one could put together a powder 



where the grains were placed in a high volume configuration, it will just sit there until 
shaken; when shaken, it will find its way to the distribution ([49|). It is possible to identify 
physical states of the powder with characteristic values of volume in our model. The 
value V = Nv± corresponds to the "deposited" powder, i.e. the powder is put into the 
most unstable condition possible, but friction holds it. When V = Nvo the powder is 
shaken into closest packing. The intermediate value of V = (vq + v\)/2 corresponds to 
the minimum density of the reversible curve. Thus we can offer an interpretation of three 

values of density presented in the experimental data [52]. Our theory gives three points 

2 



p(X = 0),p(X = oo) and p(t = 0) which are in the ratio: v , / v + v \ , v 1 1 and these 



are in reasonable agreement with experimental data: p(X = 0) = i ~ 0.64, po = ^ ~ 
0.58 &ndp(X = oo) = ^—pp^. ~ 0.62. Another important issue is the validity of the 
compactivity concept for a "fluffy", but still mechanically stable, granular array, e.g. for 
those composed of spheres with p < 0.58. In our theory, p(X = oo) corresponds to the 
beginning of the reversible branch (see Figure 11). We foresee that granular materials 
will throw much light onto glassy behaviour in the future, for questions like the pressure 
fluctuations in the jammed material are accessible in granular materials, and not entirely 
in glasses. 



5 Discussion 

In Section || we have given the analysis of stress for the case of granular arrays with 
the fixed values of the coordination number. However, a realistic granular packing can 
have a fluctuating coordination number which need not necessarily be 3 in 2-D or 4 in 
3-D. How one can extend the formalism of Section ^ to make it capable of dealing with 
arbitrary packings? The simplest (however sensible) idea is to assume that in 3-D the 
major force is transmitted only through four contacts (we call them active contacts). 
The rest of the contacts transmit only an infinitesimal stress and can be christened as 
passive. A considerable experimental evidence for this conjecture does exist. Photoelastic 
visualization experiments |5|, [2(| show that stresses in static granular media concentrate 
along certain well-defined paths. A disproportionally large amount of force is transmitted 



via these stress paths within the material. Computer simulations [31, 32, 29, 30] also 
confirm the existence of well-defined stress-bearing paths. At the macroscopic scale, the 
most obvious characteristic of a granular packing is its density. As it has been shown 
in Section |] it is natural to introduce the volume function W and compactivity X = 
^n. The statistical ensemble now includes now the set of topological configurations with 
different microscopic force patterns. Therefore, in order to have a complete set of physical 
variables, one should combine the volume probability density functional ( |47| ) with the 
stress probability density functional (|i~0[). The joint probability distribution functional 
can be written in the form 



P{V\S$] = e-f 5(V -W)H 8(sg - £ f^rf) e(f a/3 ■ n a ?) P {f a(S } . (57) 

a,P P 

We speculate that this mathematical object is necessary for the analysis of the stress distri- 
bution in granular aggregates with an arbitrary coordination number. However, this is an 
extremely difficult problem which involves a mathematical description of force chains i.e. 
the mesoscopic clusters of grains carrying disproportionally large stresses and surrounded 
by the sea of spectator particles. An explicit theoretical analysis of this problem will be a 
subject of future research. 



The aim of this paper has been to show that simple models can capture the basic 
physics of the jamming transition and provide insight into universal features of this phe- 
nomenon. We have christened the "jamming transition" as diverse physical phenomena 
which take place in granular materials, colloids and glasses. However, we argue that the 
existence of similar features like the slowing of response and the appearance of heteroge- 
neous correlated structures gives hope that there are universal physical laws that govern 
jamming in various systems. Due to the extreme complexity of this problem, the simplest 
models should be considered and solved at first. This would help to capture basic mech- 
anisms of these phenomena, and establish a theoretical framework which will incorporate 
higher complexities and details, and become a predictive tool. 
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